Quantized multi-terminal resistances in ferromagnetic graphene p-n junctions 
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The quantum Hall and longitudinal resistances in multi-terminal ferromagnetic graphene p-n 
junctions under a perpendicular magnetic field are investigated. In the Hall measurements, the 
transverse contacts are assumed to be located at p-n interface to avoid the mixing of edge states 
at the interface and the resulting quantized resistances are then topologically protected. According 
to the charge carrier type, the resistances in four-terminal p-n junction can be naturally divided 
into nine different regimes. The symmetric Hall and longitudinal resistances are observed, with lots 
of new robust quantum plateaus revealed due to the competition between spin splitting and local 
potentials. 
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Graphene, a single-layer hexagonal lattice of carbon 
atoms, has recently attracted extensive interest. 1 For 
neutral graphene, Fermi energy is located at the two in- 
equivalent Dirac points, around which graphene has a 
unique band structure with a linear dispersion, giving it 
many peculiar properties such as the half-integer quan- 
tum Hall (QH) effect. 1-4 The corresponding QH edge 
state in graphene has been proved to be chiral 2,3 with 
the chirality determined by the carrier type and the di- 
rection of the magnetic held. Thus, to manipulate the 
edge-state transmission in graphene, one may vary both 
the carrier type and concentration of graphene by the 
gate voltages 5,6 or by chemically doping the underlying 
substrate. 7 

On the other hand, a graphene p-n junction, of which 
carrier density in two adjacent regions can be individ- 
ually controlled by a pair of gate voltages above and 
below it was reported. 6 Many new and exciting phe- 
nomena are predicted, such as Klein tunneling, 8 ' 9 par- 
ticle collimation, 10 quasibound states, 11,12 and Vesclago 
lensing. 13 In addition, the electron transport through the 
graphene p-n or p-n-p junctions has stimulated many 
subsequent researches. 14-25 In the two-terminal measure- 
ments, QH plateaus show half-integer values in the case 
of unipolar junctions but fractional values for bipolar 
junctions. 6 The appearance of the fractional plateaus 
is well explained by assuming full mixing of the elec- 
tronlike and holelike edge states at p-n interface. 19 The 
experimentally observed conductance plateaus in two- 
terminal measurements are often distorted, because the 
conductance generally depends on the sample geometry 
and the interface inhomogencitics. 26-30 Thus, it makes 
sense to seek the measurements in which the QH plateaus 
are geometry-independent and do not rely on the edge- 
state equilibration at p-n interface. Different from the 
previous transport measurements, 21-25 here we propose 
multi-terminal measurements in which the transverse 
contacts are located at p-n interfaces [see Figs. 1(a) and 
1(b)]. In these proposed measurements, additional inter- 
face conduction channels appear between the transverse 
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FIG. 1: (Color online) The top panels (a) and (b) are the 
schematic for the four- and six-terminal FGPNJ devices. The 
whole devices, made of the ferromagnetic graphene, are di- 
vided into region I and region II) by the interfaces of p-n 
junctions. The bottom panel (c) shows schematically for the 
4-terminal situation the charge carrier types in both regions 
with potential difference 2V, where ni a and fi2a are the ef- 
fective chemical potentials for a spin component and Ef the 
Fermi energy. The symbol "e (h)" represents that the carriers 
are electrons (holes). 



terminals 31 and the current-carrying edge states directly 
propagate along interface without the mode mixing ef- 
fect which is unavoidable in the measurement of pre- 
vious studies. Since electron transmission between any 
two contacts is dominated by an integer number of edge 
states, the measured resistance plateaus in this situation 
should be topologically protected, which is expected to 
have a potential application in spintronic devices. 

When the Fermi energy of a ferromagnetic 
graphene 32 ' 33 is set between the spin-up and spin- 
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down Dirac points [i.e., Fermi energy lies in the spin 
gap 2M in Fig. 1(c)], the carriers consist of spin- up 
electrons and spin-down holes, each type of which is 
completely spin-polarized. 34 Therefore, one can expect 
more interesting phenomena of quantum transport in 
the ferromagnetic graphene p-n junctions (FGPNJ). In 
this paper, using the tight-binding model, we carry out 
a theoretical study on the multi-terminal resistances 
in FGPNJ in the proposed situation. Here we focus 
on the four-terminal and also discuss six-terminal Hall 
measurements, and we reveal some new fractional 
plateaus of QH and longitudinal resistances across the 
junctions. These new fractional plateaus are predicted 
to be robust against weak disorder, which is confirmed 
by a subsequent numerical calculation. 

The Hamiltonian for the four- or six-terminal FGPNJ 
in the tight-binding representation is given by, 

H = J2(^-M^- E (te^&t^+H.c), (1) 

i,<J <ij>,& 

where cj a and c^ a are the creation and annihilation oper- 
ators with spin a = ±1 at site i. gj is the energy of Dirac 
points (i.e. the on-site energy), and g; = —eV — oM 
for region I and gj = eV — <jM for region II with ±V 
the electric potentials which can be controlled by the 
gate voltages in the experiment. \i is the chemical po- 
tential. M is the ferromagnetic exchange split. 32 t is 
the nearest-neighbor hopping integral. In the presence of 
a perpendicular magnetic field B, a Berry phase factor 
4>ij = J. A' d\/4>o is added to the hopping integral with 
A the magnetic gauge potential and </>o = h/e the flux 
quantum. 

The Landau levels (LLs) of graphene in a uniform 
magnetic field B, is described by e n = ±ey(n( with 
e = VFy2eHB the energy of the first LL, vf the Fermi 
velocity, and n = 0, ±1, ±2, ... the LL index. For FGPNJ, 
the effective chemical potentials which are measured from 
the Dirac points can be expressed as /Lti CT = fi + eV + oM 
for region I and /j,2a = A f — e ^ + <J ^ f° r region II. With 
jUicr measured in unit of e, the LL filling factors can be 
given by 

v la = 2[sgn(/! K7 )(^ (T ) 2 ] + 1, (2) 

where i = 1,2, and [x] is the maximal integer no more 
than x. 

In Fig. 2, we give the energy spectrum of a two- 
terminal FGPNJ as an example to find the redistribu- 
tion behavior of filling factors owing to spin splitting and 
nonuniform local potentials. The local potential we con- 
sider here is assumed to change abruptly at the inter- 
face because the interface width can be very small. 27 ' 28 
We find that, for each spin component, there are two 
groups of LLs corresponding to region I and region II. 
By counting the number of electronlikc (holelike) LLs 
below (above) the chemical potentials (the dash lines 
in Fig. 2), it is easily to determine the LL filling fac- 
tors. In Fig. 2, for \x = 0.1, the filling factors are 
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FIG. 2: (Color online) Electron energy spectrum of a two- 
terminal FGPNJ (see Fig. 1(a) without the longitudinal ter- 
minals land 3) in the presence of a perpendicular magnetic 
field B, where the Landau gauge A = (0, Bx, 0) is adopted. 
The red and black lines denote bands of spin-up and spin- 
down states. The dashed lines represent the positions of the 
chemical potentials. Energy and chemical potentials are mea- 
sured by s = v F V2ehB with B = 10 T, M = 0.6 (~ 60 meV) 
and eV = 0.4 (~ 40 meV). 



(flt,i>i4.) = (3,-1) and {v 2 ^,v n ) = (1,-1), and for 
(i = 0.6, = (5,1) and (i^f, "24.) = (1, — 1)- 

These estimation are exactly consistent with that cal- 
culated directly by Eq. (2). Compared with the filling 
factors in pristine graphene (i/j-,1^) = (2n + l,2n+ 1) 
with n the LL index, the sequence of the filling factors 
are completely different, from which the new QH plateaus 
may be expected to arise. In the following, we first give 
the expression for the quantized resistances and then we 
make numerical calculations to confirm the predictions. 

Firstly we consider the four-terminal Hall device [see 
Fig. 1(a)]. To investigate the QH resistance #13,24 
and longitudinal resistance #13.13, which are defined by 
#13,24 = (V 2 - Vi)/h and #13,13 = (V\ - V 3 )/Ii, we 
make voltages of terminal 2 and 4 floating, and apply a 
small bias voltage Vo between terminal 1 and 3, leading 
to a steady current I\ flowing between the two termi- 
nals. From the Landaucr-Biittiker formula, 35 the electric 
current injected from terminal p is given by 

i P =jj2 T p^ v p- v iy ( 3 ) 

Here T pq is the transmission coefficient from terminal 
q to terminal p, and V q the voltage of terminal q. In 
the QH regime of the situation shown in Fig. 1(a), ex- 
cept Ti3=T3i=0, T pq between any other two terminals 
is an integer which depends on the four filling factors 
Vj^ (i=l,2, <7=t,i), since electron transmission between 
any two terminals is dominated by edge-state propaga- 
tion. Therefore, according to the positions of the effective 
chemical potentials, one can classify the system into nine 
different regimes [see Fig. 3(a)], which arc divided by 
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FIG. 3: (Color online) (a) Schematic diagram of the resis- 
tance depending on fj, and eV in FGPNJ. The resistance is 
distinguished into nine regimes which are separated by red 
solid lines and marked by different color, (b)-(g) Diagrams 
for edge states propagation in the corresponding regimes for 
V > 0, where v\ a and V2a are the filling factors in region 
I and region II. The red(black) lines with arrows represent 
spin-up(spin-down) edge states. 



four straight lines \i ia =0 (i = 1,2, <r=t,4). In FGPNJ, 
we have v\ a {^, ±V, M) = v 2<y (ii, TV, M), which is equiv- 
alent to changing p-n junction to n-p junction without 
any new physics in it. Thus, we only need to discuss the 
six regimes for V > shown in Figs. 3(b)-3(g). 

For ee-ee regime, where the first (second) "e" before "- 
" indicates the carrier type in Region I in spin-up (spin- 
down) component is electron-like while the symbol "ee" 
after "-" has a similar meaning, the nonzero transmission 
coefficients are T 2 \ = T\± = \V\\ and T 32 = T 43 = |77^"| 
with |PT| = \v\^\ + 1 1/14. |, \V 2 \ = |z/ 2 -|-| + |^24.|- The cur- 
rents I\ = —I3 = \V2\VQe 2 /h, and the voltage differences 
V x - V 3 = V and V 2 - V 4 = (\V5\/\W\)V , the QH re- 
sistance i?i3,24 and the longitudinal resistance i?i3,i3 are 
calculated (|i/i t | > \v 2 f\, \v n \ > |^|), 
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Similarly, one can obtain the resistances i?i3,24 and 
-Ri3,i3 in the other five regimes and detailed results are 
summarized as follow: 

(a) eh-eh regime (|^i t | > \u 2 f\, Wu\ < |^|): 
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13,24 



k |^lf||^2tl " 11^24-1 

e 2 h4l 2 N + MW 

h (h t | + M) 2 , 
e 2 ^24 1 2 N + Wn\W 

(b) ee-eh regime > |^2tl) : 

h \v 2 f\ 
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13,13 
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e 2 |^2il 2 + KH^I 
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Will? 



e 2 (M 2 + N|^|)N : 
(c) hh-hh regime (\u^\ < \v 21 -\, \u u \ < \v 2 i\)- 
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(e) ee-hh regime: 
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Based on the above equations, we calculate the QH and 
longitudinal resistances and plot the results for M = 
and 0.6 in Fig. 4. It displays several blocks of differ- 
ent sequence of fractional quantized resistances. In ee- 
hh regime, the QH resistances are always zero regardless 
of M, fi and V, but the longitudinal resistance are the 
summation of the contributions of region I and region II. 
Apart form the trivial cases of ee-ee, hh-hh, and ee-hh, 
the other nontrivial cases give new resistance plateaus. 
Similar to Ref. 24, both -R13.24 and i?i3,i3 have a mirror 
symmetry with respect to the /i-axis, 

-R13, 24(^14-1 v il\ ^ti ^24) = -R13, 24(^24-1 ^24; ^lfi 

-R13, 13(^14-, ^14; ^2t) ^24) = -R13, 13(^24-1 v 2i', ^lti ^14)7 

owing to the equivalence between p-n and n-p junctions, 
i.e., the equivalence between Hamiltonians H{B, )jl,—V) 
and H(B, pi, V). Furthermore, i?i3 j 24 are antisymmetric 
whereas i?i3,i3 are symmetric with respect to the V-axis, 

-^13,24(^14-! ^14 I v 2 ^,v 2 \) = 

— -Rl3,24( — ^24, — V 2 f, —^14., — ^lf), 

-R13, 13(^^14 ; ^2ti^24.) = 

-Ri3,i3(— ^24, — ^2t; — ~^it)' 
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FIG. 4: (Color online) Simulated color maps of the theoretical QH resistance -Ri3,24 [(a) and (b)] and longitudinal resistance 
-Ri3,i3 [(c) and (d)] expected from the mechanisms shown in Fig. 3 for different filling factors (^it, "14.) and (f2t; v 2i) m region 
I and region II with parameters M = and 0.6. The coordinates of the points "•" are (±M, 0) and (0,±M). (e) The one- 
dimensional slices shown in (b) (black solid) and (d) (blue solid) are plotted to show the dependence of i?i3,24 and J?i3,i3 on /1 
for a fixed eV = 0.45. 



as a result of a combination symmetry of four- 
terminal FGPNJ. This combination symmetry means 
that 1ZH(B, —fj,, y)7?. _1 is equivalent to H(B 7 fi, V), ex- 
cept the exchange of terminals 1 and 3, where 1Z = VT 
with V the particle-hole transformation \ci„ — > r)c' i(T with 
77 = 1 (— 1) if i belongs to A (B) sublattice] and T the 
time-reversal transformation. 

For M — [see Figs. 4(a) and 4(c), there is no spin 
splitting, = v\i and vi-\ = V2\\, there are only four 
different carrier-type regimes (ee-ee, hh-hh, hh-ee and 
ee-hh regimes) and the results are in exact agreement 
with recent experiment 6 and theory. 19 This indicates that 
in this situation the presence of the transverse contacts 
not only does not change the resistances, but also make 
the resistance plateaus topologically stable. For M /0, 
many new fractional values [see the values in red color in 
Figs. 4(b) and 4(d)] are induced. 39 These new fractional 
values of QH and longitudinal resistances are caused by 



spin splitting together with local potentials. More de- 
tails can be seen in Fig. 4(e), which is extracted from 
Figs. 4(b) (black curve) and 4(d) (blue curve) for fixed 
eV = 0.45. ft shows that, upon increasing the chemi- 
cal potential gradually by tuning the gate voltage, the 
QH resistance undergoes the plateau values (in the unit 
of h/e 2 ) ±1/13, ±1/9, ±1/11 and 0, where the values 
— 1/13, —1/11 and —1/9 are newly revealed, while the 
longitudinal resistance exhibits the plateau values 49/78, 
25/36, 8/11 and 1/2 with the newly induced values 49/78, 
25/36 and 8/11. We remark that since energy is mea- 
sured in unit of e, the conclusion obtained in Fig. 4 is ac- 
tually universal, i.e., it is independent of the magnitude 
of the magnetic field, so long as the system is within the 
QH regime and the magnetic length I b is much less than 
size of the device. Different from the situation discussed 
in previous studies, 6,19-23 where the edge states with dif- 
ferent filling factors are equilibrated at p-n interface, the 
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quantized resistances obtained here in the proposed situ- 
ation are topological stable, since the transverse contacts 
are just located at p-n interface and there is no mixing 
of edge states at p-n interfaces. 

Another interesting situation for four-terminal case is 
that when a small bias voltage Vb is applied to the termi- 
nals 2 and 4 instead of 1 and 3, the QH resistance i?24,i3 
and the longitudinal resistance i?24,24 can be measured. 
It is found that the QH resistances in the two cases are 
opposite to each other #13,24 = — #24,13 but the longitu- 
dinal resistances i?i3,i3 and i?24,24 are completely differ- 
ent. 

Now we explore the six-terminal Hall device [see 
Fig. 1(b)], where a small basis Vb is applied to the lon- 
gitudinal terminals 1 and 4. Here the QH and longitudi- 
nal resistances are defined by #14,26 = {V 2 — Vq)/I\ and 
#14,23 = (V2 — Vz)/I\. The resistances satisfy the proper- 
ties #14,26 = #14,35 and #14,23 = #14, 65- Similarly, there 
are nine different regimes and we still explore the resis- 
tances -Ri4,26 and #14,23 in the six regimes for V > 0, 
shown as below. A straightforward calculation gives, 

(a) ee-ee-ee regime (|z^i t | > \v 2 f\, ^14. | > |^2j,|): 
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(b) eh-eh-eh regime (\v n \ > |z^tl>Kil < I ^24.!): 



-Ki4,26 = — 1— 11.. — rrr. — i — ~. — ^ , ~ — t^a 18 -' 
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(c) ee-eh-ee regime (|^if| > l^2tl) : 

h \v 2 f\ 



(19) 



R 



14,26 — 



#, 



e 2 1^24.1(1^1 ~ |^2tl) + NN 
h 



, (20) 



V\ \ + |^24.| - |^2t| . ^ 



e 2 1^24.1(1^24.1 - ktl) + MN 
(d) hh-hh-hh regime (\u n \ < \v 2 f\, Wu\ < W21D 

R h 1 

#14,23 = 0; 
(c) eh-hh-ch regime < \v 2 ±\): 
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(f) ee-hh-ee regime: 
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FIG. 5: (Color online) The numerical QH resistance _Ri3,24 
[(a) and (b)] and longitudinal resistance i?i3,i3 [(c) and (d)] 
vs the chemical potential [i for different disorder strength W 
at M = 0.6, eV = 0.45 and 4> = 0.01. The size of the center 
block region are 9.7 nmxl0.7 nm (40 armchair chains x60 
zigzag chains). The corresponding magnetic length is Ib = 
y/h/{eB) = 0.9 nm which is much smaller than the size, so 
the edge states can propagate freely in the absence of disorder 
along the edges or interfaces without mixing with each other. 



To make sure the topological stability of the predicted 
quantized resistance, it is necessary to make a numerical 
calculation to discuss the effect of disorder on the QH and 
longitudinal resistances. 24,36 ' 37 For the four-terminal case 
discussed here, we assume that the disorder only exists 
in the center block region [dashed-box region shown in 
Fig. 1(a)]. 24,34 ' 37 In the presence of the disorder, the on- 
site energy is changed to ej + tUj , where Wi is the on-site 
disorder energy which is randomly distributed within the 
range [— W/2, W/2] with W the disorder strength. Since 
there exists no conducting channel between terminal 1 
and 3, Ti3=T3i=0. Generally, the QH and longitudinal 
resistances can be expressed as: 



#13,24 = 
#13,13 = 



T21T43 — T23X41 



T\ 2 (T 2 nTi3 + T23T4) + Ti4(T42T23 + T43T2 
T2T4 — T24T42 



T , 12(T , 24T43 + T23T4) + Ti4(T4 2 T23 + T43T2) 



(28) 
(29) 



where we introduced X4 = T41 + T42 + T43 and T2 = 
T21 + T 23 + T 24 . Here, T pq = £ CT Tr[r pff G*T ga G£] 
0, q = 1, 2, ■ ■ • ,4 and p ^ q), with T pa = z[S« - E£], 
the retarded and advanced Green's functions G^ = 
[G A ]t = [E - H c a - J2 P is the Hamiltonian 

for the center block region. The retarded self-energy 
due to the coupling to terminal p can be calculated nu- 
merically by the recutisive Green's function technique. 38 
Figures 5(a) and 5(c) show numerical results for the re- 
sistances in the absence of disorder W = 0. Compared 
with the above analytic result, these numerical results 
give exactly the same plateau values [see the plateau val- 
ues of one-dimensional slices in Fig. 4(e)]. The figures 
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also illustrate the antisymmetric and symmetric features 
of -Ri3,24 and i?i3,i3 with respect to fj, = respectively. 
In Figs. 5(b) and 5(d), we show the resistances for dif- 
ferent disorder strengths W. The plateau structures of 
■Ri3,24 and -Ri3,i3 still exist for weak and moderate dis- 
order, and the quantum plateau values are the same as 
that of the clean graphcnc p-n junction. This indicates 
that the quantum plateaus of i?i3,24 and i?i3,i3 are quite 
stable and robust against disorder due to the topological 
stability of the system. However, for a very strong dis- 
order (e.g. W > 2), the plateau structures of -Ri3,24 and 
^13,13 will finally be smeared out and then disappear. 

In summary, we have studied the QH and longitudi- 
nal resistances in the four- and six-terminal FGPNJs un- 
der a perpendicular magnetic field. In the proposed Hall 
devices, the local filling factors are independently con- 
trolled by the gate voltages and the transverse contacts 
are arranged at p-n interface. Owing to the competi- 



tion between spin splitting and local potentials, many 
new quantum plateaus of resistances are revealed. By 
avoiding the edge states mixing at the interface, the resis- 
tances are topologically protected and so robust against 
weak disorder. This multi-terminal transport study on 
the FGPNJ may provide a deeper insight into the trans- 
port properties of the FGPNJ, and the newly predicted 
fractional resistances are expected to be observed in ex- 
periment. 
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The new fractional values exclude the values in the case cussed in Ref. 34 where the resistances are i?i3,24 = ~ 

with M = and V / [see Figs. 4(a) and 4(c)], and H)/(H 2 + K| a ) and i?i 3 ,i3 = (H + H)/(H 2 + H 2 )- 

in the case with M ^ and V = which has been dis- 



